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Retrograde tracer injections in 29 of the 91 areas of the macaque 
cerebral cortex revealed 1,615 interareal pathways, a third of which 
have not previously been reported. A weight index (extrinsic frac- 
tion of labeled neurons [FLNe]) was determined for each area-to- 
area pathway. Newly found projections were weaker on average 
compared with the known projections; nevertheless, the 2 sets of 
pathways had extensively overlapping weight distributions. Repeat 
injections across individuals revealed modest FLNe variability given 
the range of FLNe values (standard deviation <1 log unit, range 5 
log units). The connectivity profile for each area conformed to a 
lognormal distribution, where a majority of projections are moderate 
or weak in strength. In the 6 2 g x 29 interareal subgraph, two-thirds 
of the connections that can exist do exist. Analysis of the smallest 
set of areas that collects links from all 91 nodes of the G 29 x 91 sub- 
graph (dominating set analysis) confirms the dense (66%) structure 
of the cortical matrix. The G 29 x 29 subgraph suggests an unexpect- 
edly high incidence of unidirectional links. The directed and 
weighted C 2 g x 91 connectivity matrix for the macaque will be valu- 
able for comparison with connectivity analyses in other species, 
including humans. It will also inform future modeling studies that 
explore the regularities of cortical networks. 
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Introduction 

A neurobiological approach to higher brain function, includ- 
ing perception, cognition, and consciousness, must be firmly 
anchored in a deep understanding of the underlying anatom- 
ical circuitry (Koch 2004; Friston 2010). Thus, it is not surpris- 
ing that understanding the connectivity of cerebral cortex 
remains a major focus of modern neuroscience (Bohland 
et al. 2009; Van Essen and Ugurbil 2012). 

Large-scale brain networks can be characterized at multiple 
levels of organization (Bressler and Menon 2010). Here, our 
level of enquiry is that of the cortical area. Localization of 
function in the cerebral cortex is closely linked to the concept 
of the cortical area and is central to our system-level under- 
standing of the neurobiology of the cerebral cortex (Rakic 
1988; Mountcastle 1997; Schiiz and Miller 2002; Rosa and 



Tweedale 2005; Zeki 2005). The pathways linking cortical 
areas have been partially elucidated by numerous anterograde 
and retrograde tracing experiments that have collectively 
revealed a number of important organizational principles 
(Zeki and Shipp 1988; Boussaoud et al. 1990; Felleman and 
Van Essen 1991; Young 1992; Goodhill et al. 1995; Barbas 
and Rempel-Clower 1997; Rockland 1997; Jouve et al. 1998; 
Sporns et al. 2000; Kaas and Collins 2001; Vezoli et al. 2004; 
Aflalo and Graziano 2011). 

Because the patterns of connections between cortical areas 
are so complex, graph theoretic approaches offer a valuable 
way to explore their network properties (Watts and Strogatz 
1998; Barabasi and Albert 1999; Watts 1999; Newman 2003; 
Boccaletti et al. 2006). Just as with many other real-world net- 
works (Boccaletti et al. 2006), cortico-cortical connections are 
directed as well as weighted. Here, we investigated both of 
these features in order to elucidate the specificity of cortical 
projection patterns. These results provide valuable infor- 
mation and constraints for future models of cortical networks. 

An early systematic study compiled the interareal connectivity 
of 32 visual areas in the macaque (Felleman and Van Essen 
1991). In the resultant 32 x 32 connectivity matrix, these authors 
reported 305 known projections out of 992 possible pathways, 
giving a link density (or graph density) of 32% for the visual 
cortex, that is, one-third of the maximum possible connections 
among visual areas were reported to exist. Subsequent studies 
using improved tracers, more systematic methods, and different 
parcellation schemes revealed numerous additional pathways. 
For example, the number of reported connections of areas VI 
and V2 expanded from a half-dozen (Felleman and Van Essen 
199D to 20 inputs to VI and 16 to V2 (Boussaoud et al. 1990, 
1991; Felleman and Van Essen 1991; Rockland and Van Hoesen 
1994; Rockland et al. 1994; Stepniewska and Kaas 1996; 
Felleman et al. 1997; Barone et al. 2000; Falchier et al. 2002; 
Clavagnier et al. 2004; Gattass et al. 2005). 

While it is well known that cortico-cortical pathways vary 
widely in their weight (connection strength), surprisingly 
little data are available on quantitative connectivity profiles in 
nonhuman primates (Falchier et al. 2002; Barbas et al. 2005; 
Burman et al. 2011; Markov et al. 2011). For example, the 
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distribution of retrogradely labeled neurons from selected in- 
jections in the macaque parietal cortex (Lewis and Van Essen 
2000) was converted into maps of labeled neuronal density 
(Van Essen et al. 2005), but not into quantitative estimates of 
interareal connection weights. Likewise, there has been little 
quantitative data on the weight consistency of any given 
pathway. A seminal paper analyzed data across studies and 
concluded that between-animal variability of a given pathway 
can exceed 2 orders of magnitude (Scannell et al. 2000). In 
contrast, a recent hemisphere-wide analysis involving retro- 
grade tracers injected into areas VI, V2, and V4 showed 
greater consistency, with interindividual variability typically 
less than an order of magnitude (Markov et al. 2011). The 
latter study also identified many newly found projections 
(NFP), especially with areas outside the classical visual 
system. Here, we apply similar sampling procedures and stat- 
istical methods in an analysis of connectivity profiles for 29 
cortical areas, thereby providing a quantitative connectivity 
database of unprecedented scope and detail. 

Previous studies reported that the vast majority of pathways 
between cortical areas are reciprocal. Felleman and Van Essen 
(1991) identified only 5 unidirectional pathways out of the 
126 explicitly tested pairs. In contrast, our analysis of a 
29 x 29 subset of cortical areas indicates an unexpectedly high 
incidence of unidirectional (nonreciprocal) pathways. 

Many existing network models assume that the network is 
sparse, that is, only a small fraction of all possible interareal 
connections actually exist (Watts and Strogatz 1998; Sporns 
and Zwi 2004; Honey et al. 2007; Hagmann et al. 2008). By 
demonstrating that the density of the binary cortical graph 
(i.e. present vs. absent) is much higher than previous estimates, 
we provide important constraints on the classes of theoretical 
networks that are relevant to modeling interareal connectivity. 

The growing efforts to explore large-scale models of cortical 
connectivity (Honey et al. 2007; Sporns et al. 2007; Rubinov 
and Sporns 2010; Adachi et al. 2012) are hindered by the fact 
that available databases have been collated from multiple 
studies using different methods of tracing, nonmatching areas, 
nomenclature and planes of section, and at best only a qualitat- 
ive assessment of connection weights (Felleman and Van 
Essen 1991; Stephan et al. 2001; Kotter 2004). Here, we report 
on a large-scale anatomical investigation of the macaque 
cortex, employing retrograde tract tracing using identical pro- 
tocols and a high-resolution analysis. This enabled us to con- 
struct an extensive and quantitative database of the weights 
and directions of interareal connections using a standardized 
parcellation scheme. Injections were made in 29 target areas 
(4 in occipital, 6 in parietal, 6 in temporal, 5 in frontal, 7 in 
prefrontal, and 1 in limbic regions). The pattern of source 
areas for each injection was determined using a parcellation of 
the entire cortical sheet into 91 architectonic areas drawn from 
published atlases (Paxinos et al. 2000; Saleem and Logothetis 
2007) and other studies. The values for the resultant connec- 
tivity profiles depend, of course, on our choice of this parcella- 
tion scheme, rather than any of the numerous alternative 
macaque cortical parcellations that are in common use (Van 
Essen et al. 2005; Van Essen, Glasser, Dierker, Harwell 2011). 
However, we believe that use of an alternative parcellation 
would have had little effect on our main conclusion (see the 
Discussion section). 

Our previous analysis of the connectivity profiles of early 
visual areas, based on statistical properties of extrinsic fraction 



of labeled neurons (FLNe) values (see the Materials and 
Methods section), included 2 important findings (Markov 
et al. 2011). 1) For repeat injections of the same area, the dis- 
tribution of FLNe values could be modeled as a negative bino- 
mial, in which the dispersion parameter determines the 
relation between the observed mean in FLN values and its var- 
iance. 2) The FLNe values of individual connectivity profiles 
followed a lognormal distribution, in which a majority of 
interareal pathways have moderate or sparse connection 
weights. Here, we demonstrate that both features, the nega- 
tive binomial and lognormal distributions, generalize to other 
cortical areas and evidently reflect important regularities of 
neocortical organization. These regularities are important 
because they provide an empirical framework for interpreting 
the results from single injections into other cortical areas. 

Our analysis of 1615 connections, including 36% that we con- 
sider to be NFP, reveals a high density of cortico-cortical con- 
nectivity (66%) and an unexpectedly high incidence (33%) of 
potentially unidirectional pathways, of which one-third (10% of 
the total) were shown to originate from corresponding locations 
within the connected cortical areas. The NFP make an impor- 
tant contribution to the connectivity profile of each area. These 
findings have important consequences for understanding corti- 
cal physiology and large-scale models of the cortex. 

The present findings provide a unique and valuable dataset 
that will aid in interpreting neuroimaging-based connectivity 
studies in humans as well as nonhuman primates. Indirect 
methods of assessing cortical connectivity using diffusion 
imaging and resting state functional connectivity are insensi- 
tive to the direction of connections and are limited in estimat- 
ing the weight of the underlying anatomical connections 
(Hagmann et al. 2008). Such data can only be obtained in 
animal models using invasive pathway-tracing methods of the 
type used here, which enable quantification of connectivity 
with single-neuron level resolution. 



Materials and Methods 

Single injections of fluorescent retrograde tracers, fast blue (FsB) and 
diamidino yellow (DY), were made in 28 macaque monkeys (27 
Macaca fascicularis and 1 Macaca mulatto). Surgical and histology 
procedures were in accordance with European requirements 86/609/ 
EEC and approved by the competent veterinary and ethical services. 
Detailed description of these methods is given elsewhere and further 
information is provided in the Supplementary material (Markov et al. 
2011). Since we use retrograde tracers, the injected area is referred to 
as the target area and the area containing labeled neurons as the 
source area. After appropriate survival times and histological proces- 
sing, high precision maps of neuron location were made using the 
software package Mercator running on ExploraNova® technology 
coupled to a fluorescent microscope stage (D-filter set 355-425 nm). 
Controlled high-frequency sampling allows stable neuron counts, 
despite the curvature of the cortex and the heterogeneity of neuron 
distribution in the projection zone (Batardiere et al. 1998; Vezoli et al. 
2004) (see Supplementary Fig. 7). Complete scanning of the hemi- 
sphere made it possible to determine the full set of ipsilateral cortical 
source areas projecting to each injected area. 

Localization of injection sites and labeled neurons was based on a 
new reference atlas that includes 91 cortical areas mapped to the left 
hemisphere of case M132. The atlas parcellation was based on a com- 
bination of histological criteria (Markov et al. 2011) and atlas-based 
landmarks (Paxinos et al. 2000; Saleem and Logothetis 2007) (see 
Supplementary material for further details). Figure 1A shows that this 
parcellation displayed on medial and lateral views of the M132 left 
hemisphere surface, generated from contours running through the 
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Figure 1. Surface atlas 3D reconstruction. {A} The 91 areas of the left hemisphere of M132 reconstructed from section drawings (see Supplementary Fig. 7). (S) Lateral and 
medial views of the 3D anatomical surface displaying areas of M132 transposed on the F99 reference brain medial and lateral view. (C) Inflated F99 right hemisphere surface, 
displaying the M132 areas registered to it. (0] F99 flat map with representation of the areas of M132. Criteria for parcellation are given in the Materials and Methods section. 



cortical midthickness (approximately layer 4). Figure IB-D shows the 
same parcellation after landmark-based surface registration to the 
macaque F99 atlas, displayed on a midthickness, inflated, and flat 
map surfaces (Van Essen 2004; Van Essen, Glasser, Dierker, Harwell 
2011). Section contours for the complete atlas parcellation are shown 
in Supplementary Figure 7. 

Quantification of Connection Weights and Identification ofNFP 

Briefly, brains were sectioned (40 [an thickness) and 1 in 3 sections 
retained for analysis. Importantly, when exploring for labeled 
neurons, 1 section in 2 (i.e. at 240 \an intervals) was examined 
throughout the cortical gray matter, which enabled the identification 
of many NFP. However, not all sections were scored when mapping 
heavily labeled areas. For details on bookkeeping of connectivity 
values, see Supplementary Materials and Methods. Areal boundaries 
on individual cases were assigned primarily using gyral and sulcal 
landmarks relative to the atlas with recourse to histological criteria 
when required (see Supplementary material). So as to derive the total 
number of labeled neurons for each pathway (see Supplementary 
Table 3), we used the observed counts generated by the frequency 
employed to estimate the expected number of neurons for 1 in 3 sec- 
tions reserved for counting. This procedure facilitates using different 
sampling frequencies for areas with very high and very low cell 



counts. Hence, the number of neurons reported in Supplementary 
Table 3 corresponds to about one-third of the expected number in the 
brain for a given injection. The uptake zones of each injection site 
correspond to a small fraction of the area injected; the volumes of the 
uptake zone are given in Supplementary Table 2. Injection sites were 
in general restricted to the cortical gray matter (see Supplementary 
Fig. 6). However, in areas VI (section 180), V4 (section 154), TEO 
(sections 187 and 195), 8m (section 161), 10 (section 20), and 7m 
(section 340), there was minor encroachment into the white matter 
(see Supplementary Fig. 6). This may have led to some contamination 
by fibers of passage, but we consider this unlikely to be a major con- 
found (see the Discussion section). 

For each injection, the number of labeled neurons in a given 
source area relative to the total number of labeled neurons in the 
brain (including those in the injected area) defines the FLN of the 
source area (Falchier et al. 2002; Markov et al. 2011). The FLNe of an 
area is estimated from the number of labeled neurons in that area rela- 
tive to the total number of labeled neurons less the neurons intrinsic 
to the injected area. For some analyses, connections are reported as 
strong (log 10 (FLNe) greater than —2), moderate (log 10 (FLNe) equal or 
less than —2 and greater or equal to —4), or sparse (log 10 (FLNe) less 
than -4). 

To identify connections to the 29 target areas that were not pre- 
viously reported, we assessed the connections reported in 119 
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published studies of connectivity (see Supplementary Tables 1 and 6). 
We adopted a conservative position with regard to categorizing the 
NFP. In some publications, connections that we report as known were 
illustrated but not explicitly listed in the text of the publication. Many 
such connections are not represented in existing databases. Here, 
only projections that previously were neither depicted nor listed in 
publications or databases (cf. see Supplementary Tables 1 and 6) are 
classified as NFP. The veracity of labeled neurons were checked by 
senior team members (H.K., CD., P.M., N.T.M., and C.L.). 



Modeling Inconsistency 

It is critical to evaluate the reliability of observing a given connection 
across multiple injections. We refer to this as the consistency of a con- 
nection or pattern of connections. Repeated injections into areas VI, 
V2, V4, and 10 were used to assess the consistency of projections to 
the early visual areas and a representative prefrontal area. In our pre- 
vious study, the negative binomial model provided a good description 
of the variability of the mean counts of the inputs to the early visual 
areas (Markov et al. 2011). Using identical methods to evaluate the 
data from the area 10 injections, we confirmed that area 10 is also 
well described by the same negative binominal model. By extension, 
we assume that a negative binomial model with similar dispersion 
should be valid for injections in other regions. For a negative bino- 
mial distribution with a known dispersion, there is a functional 
relation between the variance and the mean, that is, the variance of 
the negative binomial distribution is ft + ft 2 / 0, where pi is the mean 
and 9 the dispersion. We used this relation in order to estimate the 
precision of projections from individual injections as well as to model 
the consistency of the sparse connections. 

To analyze statistical characteristics of inconsistent projections, we 
calculated the probabilities of observing no neurons under several 
models. For the simple case of the Poisson distributed counts, the 
probability of observing y counts given the mean number of counts, 
fj, is: 



For y — 0, 

P(y = 0| M ) = (1) 
For the negative binomial distribution, the probability density is: 



P(y\^ 0) = 



T(y + v) ( 6 YY M V '' 



r(e)ro+ 1) \o+ti. 



where T is the gamma function and 6 the dispersion. Then, the prob- 
ability of observing zero counts is 



P(y = 0| M , 0) = 



(2) 



For the special case of the geometric distribution (d = 1), the prob- 
ability of zero counts is: 



P(y = 0\pi, ff) = 



1 + M 



(3) 



Let p be the probability of observing zero counts from a projection 
from a single injection, as given by any of equations (l)-(3). The 
probability of observing at least one or more neurons in a projection 
(i.e. y>0) from a single injection is 1 —p. Since each injection is inde- 
pendent, the probability that some neurons are observed in each of n 
injections, is (1 —p}", and the probability of observing no neurons in 
at least 1 of n injections is 1 — (1 —/>)". This represents the probability 
of observing zero counts in one or more of the n injection exper- 
iments performed. 

Updates, atlases, and additional information are available at 
www.core-nets.org. Surface-based atlas datasets are accessible at 



http://sumsdb.wustl.edu/sums/directory.do?id=8287442&dir_name= 
MARKOV_CC12. 



Results 

Injections were included in our analysis if they showed little 
or no involvement of the underlying white matter and were 
restricted to a single cortical area. The extent of each injection 
relative to areal borders and white matter is shown in line 
drawings of selected section contours for each individual case 
in Supplementary Figure 6. The full set of 39 injections in 29 
areas after mapping to the atlas surface is shown in Sup- 
plementary Figure 1. This mapping provides each injection 
site with a well-defined stereotaxic location that is indepen- 
dent of the underlying cortical parcellation for the stereo- 
taxic coordinates of the injection sites (see Supplementary 
Table 7). Hence, while the focus of our analysis is on the 
pattern of area-to-area connectivity, these data can also be 
evaluated in terms of the areal inputs to a set of atlas coordi- 
nates that have well-defined locations relative to alternative 
parcellation schemes that have been mapped to the F99 atlas 
(Van Essen, Glasser, Dierker, Harwell 2011; see the Discus- 
sion section). 

We first show how modeling the observed variability of 
FLNe in repeat injections allows the estimation of the 
reliability of results obtained from a single injection. One 
important question is how closely the weight of a pathway 
observed after a single injection approximates the average 
connection weight determined from repeat injections. A 
second issue concerns the consistency of connections. If a 
given interareal pathway is observed after a single injection, 
what is the likelihood that the same pathway will be observed 
after repeat injections, and to what extent does this likelihood 
depend on the observed connection weight (FLNe)? To 
address these issues, we extended our previous analysis of 
connectivity variability for visual areas VI, V2, and V4 
(Markov et al. 2011) to include prefrontal area 10, which (as 
shown below) is more extensively connected than are early 
visual areas and might in principle show different variability 
dependencies. 

Modeling Variability of FLNe Values in Area 10 

Repeat injections were aimed at area 10 at the very rostral tip 
of the cortex. Figure 2 shows four 3D views of each recon- 
structed injection site in relation to the pial surface (gray) and 
the gray/white border (blue). Inspection of these reconstruc- 
tions, especially the medial views (far right column), confirms 
that the uptake zones were restricted to the frontal pole, span- 
ning the cortical gray matter but not encroaching on the 
underlying white matter (see Markov et al. 2011 for details 
on uptake zone determination). The 2 smallest injections 
(M131LH and BB341LH, Fig. 2A,D) had almost identical 
locations in the rostral pole. The 2 larger injections (M136LH 
and M137LH, Fig. 2B,C) also included the anterior pole but 
extended further ventro-laterally. 

The 4 repeat injections in area 10 showed similar labeling 
connectivity profiles in terms of the identity of labeled areas 
and the FLNe values for each projection. In Figure 3^4, the log 
(FLNe) values are ordered by the geometric means of neuron 
counts (thin red line) for all areas that project to area 10 in 
one or more cases. Fifteen projections had a mean FLNe value 
exceeding 10~ 2 , which we classify as a strong projection. 
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Figure 2. Three-dimensional reconstruction of the 4 injection sites in area 10 (respectively, cases 34, 37, 38, and 39). Injection site in red and white matter in blue. Dotted line 
shows the limits of area 1 0. (i) Fronto-lateral view, (ii) frontal view, (iii) fronto-medial view, and (iv) medial view. 



Twenty-seven projections had FLNe values between 10~ 2 and 
10~ 4 , which we classify as moderate in strength. Twenty-four 
projections had FLNe values below 10~ 4 , which we classify as 
sparse. The geometric mean is reasonably well fit by a lognor- 
mal distribution (thick black line), as previously observed for 
visual areas VI, V2, and V4 (Markov et al. 2011). 

The observed range of FLNe values for the 4 repeat injec- 
tions in area 10 is less than one order of magnitude for most 
areas. However, for 4 input areas (45 A, 31, 7 A, and DP) the 
range exceeds 2 orders of magnitude. In Figure 3B, the stan- 
dard deviation (SD) and mean is fit by a negative binomial 



distribution (green line) having a dispersion parameter of 5.0 
(see the Methods section). The best-fitting negative binomial 
curve differs considerably from a Poisson distribution (red 
line) except for low FLNe values. It is closer to a geometric 
distribution (blue line), as was seen for injections of areas VI, 
V2, and V4 (Markov et al. 2011). The dispersion parameter 
obtained from area 10 (5.0) yields a curve similar to that 
obtained from the visual areas (average dispersion, 7.6). 

Some of the observed variability might be attributable to 
measurement errors owing to imperfect delineation of areal 
boundaries for the various source areas. One way to assess 
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Figure 3. Variability of labeling after repeat injections in area 10. Ifi) FLNe plotted by area for 4 injections in cortical area 10 ordered by the geometric mean (thin red curve) of 
the values for each projection (excluding data points for "absent" projections when computing the mean). The thick black curve indicates the expected values for an ordered 
sample from a lognormal distribution with the same mean and SD. Symbols below 10~ 6 indicate zero values. (3-D) The SD as a function of the mean. The curves are the 
predictions for a Poisson (red), geometric (blue), and the best-fitting negative binomial distribution (green). The dispersion parameter of the negative binomial distribution and its 
95% confidence interval are indicated in the inset; (fl) values for the 4 injections in area 10; (C) areal values for the repeat injections in V1, V2, V4, and area 10; (D) cumulative 
regional values. 



variability due to errors in the location of area borders is to 
compare the variability of the repeat injections in areas VI, 
V2, V4, and 10 (Fig. 30 with the variability observed when 
the areas are grouped into much larger regions (Fig. 3D). 
Comparisons of the slopes in Figure 3C,D show only a 
modest reduction in variability, suggesting that areal delinea- 
tion is not the major source of variability. 

In terms of the presence versus absence of connections, 
the 4 repeat injections in area 10 showed similar labeling pat- 
terns (Fig. 3A). For the 67 areas containing label from at least 



one injection, 50 (75%) were labeled by all 4 injections, 7 by 3 
injections, 4 by 2 injections and 6 by only a single injection; 
24 areas lacked connections to the injection area in all 4 
cases. Of the 17 areas labeled by fewer than 4 injections, 15 
were very sparse, with labeling of fewer than 10 neurons on 
average. These were near the limit of systematic detectability 
expected for a variable where the heteroscedasticity of the 
SD/mean relation follows a negative binomial distribution 
(Fig. 3B). The 2 exceptions were areas 7 A and DP, with mod- 
erate labeling following a single injection. 
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Figure 4. Relationships of means and 95% confidence intervals from multiple injections to the values from single injections. (4) Histogram of residuals for multiple injections 
(VI, V2, V4, and 10) with respect to lognormal order statistics normalized to unit area. Dashed curve: A kernel density estimate of the underlying distribution obtained by 
convolution of the histogram with a Gaussian; solid curve: Best-fitting normal distribution (mean = -0.003, SD = 0.430). (6) Ordered FLNe values from a single injection in 
cortical area V1 (white circles) with 95% confidence intervals expected on the basis of a negative binomial distribution (error bars). The small black dots correspond to values 
obtained from 4 other injections in the same area. The blue circles are the geometric means. For the 3 entries on the far right (MB, 8r, and 7op), there were no labeled neurons 
from the V1 injection used for FLNe rank ordering. (C) Ordered FLNe values from a single injection in cortical area 1 0 (white circles) with confidence intervals and small black 
dots (3 other injections) and blue circles as described for area V1 . For the 8 entries on the far right, there were no labeled neurons from the area 1 0 injection used for FLNe rank 
ordering. 



Using Single-Injection FLNe Values to Predict the Mean 
from Multiple Injections 

For the multiple injections in VI, V2, V4, and area 10, the 
histogram in Figure 44 shows the residual between the log 
(FLNe) values from individual injections and the predicted log- 
normal curve. This provides a quantitative measure of the 
relationship between average connection weights and the 
values obtained from individual injections. Given an FLNe 
value from only a single injection, the mean value (i.e. the value 
for the lognormal distribution) lies within a factor of 2.7 (1 SD) 
in about 68% of cases. The observed variability presumably 
reflects a combination of factors (see the Discussion section). 

The analysis illustrated in Figure 4B,C and Supplemen- 
tary Figure 3 indicates how accurately FLNe values from a 
single injection can predict the mean obtained from multiple 
injections. Figure 4B shows the ordered FLNe values from a 
single VI injection (white circles), with a 95% confidence 



interval assuming a negative binomial distribution analysis 
with the dispersion observed following the multiple injec- 
tions. The difference between the observed single-case value 
and the mean of 5 injections (blue circles) is generally quite 
small and for this exemplar injection in every case lies within 
the 95% confidence interval predicted by the negative bino- 
mial distribution. The maximum difference is 0.99 log unit, 
corresponding to a 9-5-fold difference between the single- 
subject and group-average results. A similar analysis for area 
10 (Fig. 4C) shows a comparably good fit for most of the data 
involving strong and moderate projections. However, for a 
few projections, the difference between individual and group 
average exceeds an order of magnitude. For the full set of 14 
repeats (see Supplementary Fig. 3), involving 544 single FLNe 
values, the 95% confidence values included 98.5% of the 
observed FLNe means determined from repeat injections. 
Altogether, this analysis indicates that FLNe values obtained 



Cerebral Cortex January 2014, V 24 N 1 23 



from single injections are usually within a factor of 3 and 
highly likely to be within a factor of 10 of the mean value. 
However, exceptions can occur, especially for sparse connec- 
tions, as discussed in the next section. 

Consistency of Sparse Connections 

Elsewhere we have shown that increasing the sampling rate 
reduces the variability in estimates of the FLNe (Vezoli et al. 
2004). Here, we consider how sampling rates influence the 
consistency. 

For the 37 projections having interpolated mean counts 
fewer than 10 neurons, 31 were inconsistent (see Supplemen- 
tary Table 3). Here, we want to distinguish between inconsis- 
tency attributable to sampling error (e.g. due to incomplete 
sampling given that not all sections were examined) and 
genuine biological variability in which some pathways are 
present in some but not all cases. The following analysis indi- 
cates that statistical fluctuation due to sampling errors can 
largely account for the observed inconsistency of very sparse 
projections. 

Figure 5A shows the probability of observing zero counts 
(i.e. failing to detect a connection that exists, namely a false- 
negative) as a function of the mean (expected) number of 
labeled neurons, f,i, after a single tracer injection for Poisson, 
negative binomial, and geometric distributions (eqs 1, 2, and 
3, respectively, in the Materials and Methods section). The 
probability of a false-negative falls below 0.05 (grey horizon- 
tal line in Fig. 5A) for means as low as 4 neurons for both the 
negative binomial and Poisson curves, consistent with the 
convergence of their variance/mean ratio at small FLNe 
values shown in Figure 3B. The probability of observing at 
least one false-negative decreases is a sigmoid function of the 
mean, as shown in Figure 5B for 2, 3, and 5 injections (span- 
ning the range in our data set). For 5 repeat injections, as in 
our VI injections, the false-negative probability falls below 
0.05 (grey horizontal line) for a mean of 7 neurons. The prob- 
ability of observing at least one false-negative after 5 injec- 
tions drops much more steeply for a negative binomial than 
for a geometric distribution (Fig. 5C7). 

Taking the observed mean as the best estimate of the popu- 
lation mean (or expected value), the probability of the 
observed inconsistency under the negative binomial model 



exceeds 0.05 for 27 of 39 projections having an observed 
mean fewer than 10 neurons (see Supplementary Table 3, 
column G). Of the 12 cases with probabilities below 0.05, 3 
are consistent and had an observed mean exceeding 7 
neurons. After applying the Bonferroni correction for evaluat- 
ing multiple probabilities (Bretz et al. 2010), only the projec- 
tions of areas 7A and DP to area 10 were inconsistent despite 
a statistically very low probability of being so. Genuine bio- 
logical variability presumably exists in the connectivity pro- 
files for different individuals, and it is possible that some 
pathways are present in some individuals and altogether 
absent in others. However, statistical fluctuations in the data 
largely accounts for the observed inconsistency of very sparse 
projections (see the Discussion section). These observations 
on inconsistency depend on the fine sampling that we em- 
ployed. In simulations of sparse projections having a negative 
binomial distribution, coarser sampling had little effect on the 
expected proportion of zero counts observed, but the variabil- 
ity in the estimates was proportional to the square root of the 
sampling interval (e.g. a 1:16 sampling interval results in 4 
times as much variability in the estimated occurrence of zero 
counts compared with a 1:1 sampling). 

In summary, this analysis of repeat injections provides ob- 
jective constraints on what can and cannot be learned from 
making a single injection into any given cortical area. Specifi- 
cally, a single-injection approach can 1) detect all but the 
most sparse projections with high probability, 2) provide a 
reasonable estimate of the connection weight of each pathway 
(generally within an order of magnitude), and 3) identify 
some of the sparse connections that are statistically likely to 
be inconsistent across multiple injections. Repeat injections 
enable identification of a few additional sparse projections 
and also provide better estimates of average connection 
weights. However, given the paucity of quantitative data on 
interareal connection weights in the macaque (see the Intro- 
duction section), the single-injection results described below 
have advantages in terms of the overall information gained 
when applied to a large number of areas. 

Connectivity Patterns Revealed by Single Injections 

Individual injections made into 25 cortical areas (in addition 
to areas VI, V2, V4, and 10) provided extensive new 
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Figure 5. Theoretical analysis of projection consistency. (4) Probability of observing zero counts as a function of the true mean for Poisson (red), negative binomial (green), and 
geometric (blue) distributions. (6) Probability of observing at least one case of zero counts in n injections as a function of the true mean for a negative binomial distribution with 
dispersion parameter equal to 7.2. (C) Comparison of the probability of observing at least 1 zero as a function of the true mean in n = 5 replications for the geometric (blue) and 
negative binomial [6 = 7.2) distributions. 
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information of several different types. This includes the identi- 
fication of many NFP, quantification of the weight of all path- 
ways (based on FLNe values), demonstration of the lognormal 
distribution of FLNe values as a general principal of the 
cortex, and evidence for many unexpectedly unidirectional 
rather than reciprocal pathways. 

We begin by discussing results for an exemplar injection of 
area F2 in the motor cortex. Figure 6 shows selected section 
contours with the injection site in solid red, the injected area 
(F2) in gray, and retrogradely labeled neurons as red dots; 
labeled neurons intrinsic to the area are not shown. Areal 
boundaries are indicated by black lines. Source areas repre- 
senting NFP are labeled in bold blue and marked by an aster- 
isk; the sections were selected to illustrate all the NFP. In a 
few places, the label in a restricted region (identified by a 
black rectangle) comes from an adjacent section projected 
onto the section shown. Red lettering indicates inferred 
known projections, which although illustrated in previous 
studies listed in Supplementary Table 6 were not explicitly 
reported largely due to differences in parcellation scheme 
used (see below). 

Area F2 is located in the dorso-caudal part of Brodmann's 
premotor area, area 6. It is located anterior to Fl (area 4), 
extends rostrally up to 3 mm in front of the genu of the 
arcuate sulcus and is bordered ventrally by the spur of the 
arcuate sulcus (Matelli et al. 1998). In this case, the injection 
was dorsal and rostral to the superior precentral dimple 
(Fig. 6 and see Supplementary Fig. 7) and was in the hind 
limb representation (Dum and Strick 1991; He et al. 1993; 
Godschalk et al. 1995; Graziano and Aflalo 2007). 

Fifty one projections were identified from this injection, of 
which 14 we consider NFP. The known projections include 6 
areas in the limbic cortex (areas 23, 24a, 24b, 24c, 24d, and 
31) (Barbas and Pandya 1987; Ghosh and Gattera 1995; 
Matelli et al. 1998; Caminiti et al. 1999; Petrides and Pandya 
1999; Marconi et al. 2001; Tanne-Gariepy et al. 2002; Luppino 
et al. 2003; Morecraft et al. 2012), 8 areas in the prefrontal 
cortex (areas 8B, 81, 8m, 8r, 44, 46v, 9/46d, and 9/46v) 
(Barbas and Pandya 1987, 1989; Ghosh and Gattera 1995; Ca- 
miniti et al. 1999; Petrides and Pandya 1999; Marconi et al. 
2001; Luppino et al. 2003; Takada et al. 2004; Petrides and 
Pandya 2006; Morecraft et al. 2012), 8 areas in the frontal 
cortex (areas Fl, F3, F4, F5, F6, F7, INSULA, and parainsula) 
(Barbas and Pandya 1987; Ghosh and Gattera 1995; Caminiti 
et al. 1999; Petrides and Pandya 1999; Marconi et al. 2001; 
Tanne-Gariepy et al. 2002; Luppino et al. 2003; Takada et al. 
2004; Morecraft et al. 2012), 11 areas in the parietal cortex 
(areas 3, 5, 7A, 7B, 7m, 7op, AIP, LIP, MIP, VIP, and SII) 
(Jones et al. 1978; Petrides and Pandya 1984; Ghosh and 
Gattera 1995; Johnson and Ferraina 1996; Matelli et al. 1998; 
Caminiti et al. 1999; Cipolloni and Pandya 1999; Petrides and 
Pandya 1999; Luppino et al. 2001; Marconi et al. 2001; Tanne- 
Gariepy et al. 2002; Petrides and Pandya 2009; Morecraft et al. 
2012), and 4 areas in the temporal cortex (areas MST, STPc, 
STPi, and TPt) (Seltzer and Pandya 1989; Luppino et al. 
2001). Previous evidence for 4 of these projections is rela- 
tively sparse: TPt (Luppino et al. 2001), 7op (Ghosh and 
Gattera 1995; Caminiti et al. 1999; Cipolloni and Pandya 
1999), and 9/46v (Takada et al. 2004), so strictly speaking, 
they could also be considered inferred rather than known. 

NFP are located in temporal (9), prefrontal (4), and frontal 
(1) regions (blue areal labels in Fig. 6). These include area MT 



(Fig. 6A), area MB (Fig. 6D), area PERHIRHINAL (Fig. 6E), 
area PBr (Fig. 6F), areas IPa, TEa/ma (Fig. 6G), area PGa 
(Fig. 6H), area STPr (Fig. 61), area TEMPORAL_POLE 
(Fig. 6J), area ProM (Fig. 6X), area OPRO (Fig. 6Z), area 45B 
(Fig. 6M), and areas 12 and 13 (Fig. 6N). These NFP show 
relatively few labeled neurons in the sections illustrated. 
Quantification of each pathway (see Supplementary Table 6) 
indicates that most NFP to F2 are indeed sparse, but several 
(areas 12, MB, and PGa) are moderate in strength. 

Different subregions of F2 have distinct patterns of inputs 
(Johnson et al. 1996; Matelli et al. 1998; Tanne-Gariepy et al. 
2002; Luppino et al. 2003), so it is important to specify the 
location of the injection site. For example, the injection in 
case 27 was in the hindlimb representation and did not lead 
to labeling in areas V6A, 32, and PIP, whereas more lateral 
and ventral injections (in the arm representation) do label 
these areas (Johnson et al. 1996; Matelli et al. 1998; Caminiti 
et al. 1999; Petrides and Pandya 1999; Tanne-Gariepy et al. 
2002). 

Figure 7 shows a connectivity map for the F2 exemplar in- 
jection displayed on cortical surface maps (lateral and medial 
inflated maps plus a flat map of the F99 atlas). Connection 
strengths are encoded as sparse, moderate, or strong (dark to 
light shades) using green for previously reported projections 
and red for NFP. Several general observations emerge from 
examination of the F2 connectivity map (Fig. 7) and the other 

28 maps (data not shown). In general, the areas providing 
inputs to any given target area form a group that covers at 
least half of the cortical surface (e.g. area 24c) and in some 
cases more than 3 quarters of the hemisphere (e.g. area 7A). 
The strongest connections tend to be from areas close to the 
injection, but there are exceptions. NFP are rare in the 
immediate vicinity of the injected area, but are common at 
both intermediate and distant locations and often form a near- 
continuous belt of areas, as in Figure 7. 

Weight Distribution of Known Projections and NFP 

Injections in all 29 target areas revealed a total of 1,615 corti- 
cal pathways. Of these, 579 have not been previously de- 
scribed (NFP). Each target area received projections from 
between 26 to 87 source areas. Supplementary Figure 2 
shows the precise location of injection sites and the pattern of 
projections in selected section contours for 28 target areas (all 
but F2). Supplementary Table 6 provides citations for each of 
the 1,036 known projections reported in this study for all the 

29 target areas. The evidence for the existence of the known 
connections is generally relatively good, insofar as most pro- 
jections are reported in several studies. The incidence of in- 
ferred projections (suggestive but not compelling evidence) 
described above for area F2 was similar for the other 28 
target areas. 

For all 29 areas injected, including the 14 repeat injections 
in 4 areas, we found that the ordered FLNe was well described 
by a lognormal distribution (Fig. 8 and see Supplementary 
Fig. 4). The curves in each plot are based only on the mean 
and SDs of the loglO (FLNe) values and the number of projec- 
tions; no free parameters were used to constrain the shape of 
the curves. This consistency in the weight distribution indi- 
cates a strong regularity in the cortex. The profiles shown in 
Figure 8 and Supplementary Figure 4 include 95% confidence 
intervals based on the negative binomial model fit to the data 
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Figure 6. Charts of labeled neurons following injection in area F2. Upper left: Section levels ffl-O) indicated on a lateral view of the cortex, red filled region indicates pick-up 
zone of injection site. (4-0) Charts of coronal sections of retrogradely labeled neurons (red dots). Black rectangle indicates neurons from nearby sections. Blue lettering and 
asterisk identifies NFP Red lettering identifies inferred known projections (see text). Scale bar: 2 mm. 



and constitute the connectivity profiles of the 29 injected 
areas. These curves are necessarily monotonically decreasing 
by the way they were constructed. However, there is no 
a priori reason why the distribution must be well fit by a 



lognormal distribution, and we consider this an important 
experimental finding. The curves shown in Figure 8 illustrate 
common features in the connectivity profiles including just 
how closely the single FLNe values track the lognormal 



26 Density of the Cortical Network 



Markov et al. 




Sparse Medium Strong 



Known 
NFP H M 

Figure 7. Cortical surface maps for the F2 exemplar injection. Flat map plus medial and lateral inflated maps for injections shown in Figure 8. Connection strengths are encoded 
as sparse, moderate, or strong (dark to light shades) using green shading for previously reported projections and red for NFP The area injected is in black. 



distribution as well as the progressive widening of the 95% 
confidence interval with decreasing FLNe values. 

Supplementary Table 1 lists the NFP together with the 
known pathways and supporting publications for each target 
area. For the target areas in each cortical region, the NFP con- 
stituted a substantial fraction of the total number of connec- 
tions (limbic 13%, prefrontal 40%, frontal 30%, parietal 30%, 
temporal 43%, and occipital 37%). Figure 9 shows a histogram 
of connection weights for known projections (white bars) and 
NFP (red bars) in intervals of 0.5 log 10 (FLNe). Known connec- 
tions are on average stronger than NFP, but the 2 populations 
overlap extensively. Remarkably, 43% of the NFP had FLNe of 
moderate strength, and a few (2%) are classified as strong con- 
nections. For very low FLNe values, NFP correspond up to 
90% of the population, but constitute a decreasing fraction of 
the source areas with increasing FLNe. 

Multiway contingency tables make it possible to compare 
the consistency of NFP and known projections (see Sup- 
plementary Fig. 5). This shows that at similar weights, the con- 
sistency of the NFP is similar to that of the known projections. 

Reciprocity 

Connections between any pair of areas can be categorized as 
bidirectional (reciprocal), unidirectional (a connection ob- 
served in only one direction), or unconnected. Previous ana- 
lyses of this issue have been hampered by the incompleteness 
of data testing for connectivity in both directions using con- 
sistent parcellation criteria (see the Introduction section). The 
availability of 29 injections analyzed using the same parcella- 
tion scheme enabled us to test for the existence of all path- 
ways and to generate an edge-complete graph, which we call 
G29x29- In our (£29x29 edge-complete matrix, the connectivity 
was bidirectional between 214 projection pairs, unidirectional 
for 108 pairs, and unconnected for 84 pairs (see Supplemen- 
tary Table 4). This is a much higher incidence of apparent 
unidirectional connections than reported in previous analyses 
(see the Discussion section). Of the 108 unidirectional con- 
nections, 4 are classified as strong and 36 as medium (see the 
Materials and Methods section). For bidirectional connections, 
we compared the ratio of the FLNe weights for each direction 
in order to assess the degree of symmetry (see the Materials 
and Methods section). We categorized 126 reciprocal path- 
ways as symmetric and 88 as asymmetric (see Supplementary 
Table 4). 

As already noted, some cortical areas have nonuniform 
connectivity across their extent, and this could impact assess- 
ments of reciprocity. For example, whereas foveal VI and V4 



are reciprocally connected (Zeki 1978; Zeki 1980; Van Essen 
et al. 1986), peripheral V4 projects to VI but does not receive 
a projection from VI (Kennedy and Bullier 1985; Perkel et al. 
1986). We therefore investigated whether our observation of 
unidirectional connections may in some cases reflect a misa- 
lignment of local connectivity patterns for pathways that are 
actually bidirectional in an area-to-area sense. Consider a 
putative unidirectional projection from area X to area Y. We 
examined the location of labeled neurons in area X (resulting 
from the injection in area Y) with respect to the injection site 
in area X (from a different case) that failed to label neurons in 
area Y. If the injection site of area X overlaps with the labeled 
neurons projecting to area Y, we consider this strong evidence 
for a bona fide unidirectional pathway. For example, 
Figure 10^4 shows labeled neurons in area 9/46v following in- 
jection in TEO (case 11), and in gray the injection site location 
in area 9/46v (case 31) that failed to label neurons in TEO. 
Figure 10B-E shows 4 additional examples of overlap 
between unidirectional projections and the injection site 
testing the reverse direction. For the 58 out of 108 apparently 
unidirectional pathways in which more than 10 neurons were 
labeled, 32 showed an overlap between the injection site and 
labeled neurons indicative of a genuine unidirectional projec- 
tion (see Supplementary Table 4). This analysis suggests that 
a minimum of 10% of the 322 connected pairs are genuinely 
unidirectional, and over half of these have medium to strong 
FLNe values (see also the Discussion section). 

A 29 x 91 Weighted Connectivity Matrix 

Figure 11 A shows the 29x91 connectivity matrix in which 
the color of each entry represents the loglO (FLNe) value for 
that pathway (brighter shades representing stronger connec- 
tion weights; black represents no connection). Each column 
gives the FLNe profile of inputs observed for a given area and 
each row its outputs. The matrix is asymmetric, in contrast to 
many connectivity matrices obtained using neuroimaging 
methods that cannot distinguish the directionality of connec- 
tions. The rows and columns were ordered so as to maximize 
the overall similarity between neighbors (see the Materials 
and Methods section). Visual areas are concentrated in the 
upper left quadrant of the matrix. Motor and somatosensory 
areas are concentrated in the lower right quadrant. Higher- 
level areas of the frontal, parietal, and temporal cortex are 
mainly in the middle portion. Green squares indicate the 
"identity" entry for the same area in a row and column. These 
entries run approximately along the diagonal, suggesting that 
areas having a similar pattern of inputs also have a similar 
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Figure 8. Connectivity profiles for 6 injected areas, chosen to illustrate a greater than 3-fold range in number of connections (in-degree distribution, see Fig. 12). The log(FLNe) 
values are ordered. The solid curves correspond to the predicted order statistics for a lognormal distribution with the same mean and SD as the data. The error bars are 95% 
confidence intervals, assuming that the data follow a negative binomial distribution with dispersion equal to 7. Connectivity profiles for the remaining injections are displayed in 
Supplementary Figure 4. 



pattern of outputs. However, the many deviations from the 
diagonal suggest that some areas having similar input profiles 
have distinctly different output profiles, and vice versa. 
Figure 11B shows a 29x29 matrix that includes only source 
areas that also served as injection sites. Deviations from sym- 
metry along the diagonal for black entries signify asymmetries 
in the existence of observed projections. Deviations from sym- 
metry in the shaded entries signify asymmetries in the 
strength of these projections. 

Graph Density 

Existing databases do not provide reliable estimates of the 
density of the network of cortical areas, nor has there been a 



concerted effort to extrapolate from the existing data to the 
connectivity of the full interareal network (FIN). Here, we 
explore the classical notion of graph density (Janson et al. 
2000; Newman 2010) of the FIN using multiple approaches 
that give converging results. This analysis draws from the 
weighted connectivity matrix (Fig. 11) using binary measures 
of connectivity (i.e. connections existing or not, independently 
of their strength). Graph density is a fundamental measure of 
the graph's overall connectedness, extensively used in 
network science and also in earlier analyses of cortical connec- 
tivity (Sporns and Zwi 2004; Bullmore and Sporns 2009). 

Based on the M132 atlas parcellation (Fig. 1), the FIN con- 
tains 7V FIN = 91 cortical areas that represent the nodes of the 
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Figure 9. Weight comparisons for known projections and NFP Distribution of known 
projections and NFP as a function of projection magnitude (FLNe) at intervals of 0.5 
log 10 , following the injection of the 29 target areas. Blue line indicates the 
percentage of NFP within each interval. 



G 91x91 graph. The directed edges of the FIN correspond to 
directed connections between nodes, based on the FLN. Our 
analysis of the FIN makes use of the G29X91 directed sub- 
graph of projections within FIN, which reveals all the in- 
degrees of the injected 29 nodes. It also makes use of the 
G29 x 29 edge-complete subgraph of FIN, corresponding to the 
connections among just the 29 injected areas. Both G 2 9*9i 
and G29 x 29 contain complete information about the status of 
their edges and would not be influenced by injections into 
additional areas elsewhere in the cortex. Given that the 29 
injected areas are widely distributed among the 6 regions, 
the G 29 x 29 subgraph is likely to reflect major characteristics of 
the FIN. 

The density of a directed graph is given by the ratio p = M/ 
[N{N—Y>] between the number of directed edges (links) M of 
the graph and the total number of possible links, N{N—1), 
where N is the number of nodes in the graph. The G29 x 29 
graph has M= 536 (binary) directed links from the maximum 
possible of N(N— 1) = 812, and therefore, it is strongly inter- 
connected, with a graph density of p = 0.66 (66%). Because it 
is an edge-complete subgraph of FIN, the density of G29 x 29 is 
expected to be comparable to that of the FIN. 

The in-degrees of the G 2 9x9i graph (i.e. the number of 
source areas projecting to each of 29 target areas, Fig. 12) range 
from 26 to 87 with a mean of (k) in = 55. 4; their distribution 
(right side of figure) is concentrated around the mean. The 
density of the FIN was estimated as follows. Because every 
directed edge is an in-link to some node, the total number of 
edges M FIN equals the total number of in-links in the FIN. We 
lack data on the in-links to nodes that were not injected, but we 
can assume that they are characterized by the same average 
in-degree as the 29 injected nodes. Assuming M FIN ~ (k) in 
■Wfin = 5,071 for the FIN ((£91x91) leads to the prediction 
p FIN = M FIN /[Af FIN (7V FIN -l)] w (k)J(N FlN -l) « 0.62 (62%), 
which is of the same order as the density for the edge- 
complete graph G 29 x 29 . 

A dominating set analysis on G 29x91 provides further evi- 
dence that the FIN is indeed dense. In graph theory, a subset 
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Figure 10. Positive identification of unidirectional pathways. Examples of contiguity 
of injection site and retrograde labeled neurons in unidirectional pathways. For a 
unidirectional projection X>Y we show the labeled neurons in area X following 
injection in area Y. The injection site in area X is the injection site that failed to label 
neurons in area Y. Note that the injection site and the reported labeled neurons are in 
2 different brains. Scale bars: 2 mm. 
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Figure 11. Weighted connectivity matrix. (4) Each row represents 1 of the 91 source areas; each column represents 1 of the 29 injected target areas. The color shows the 
strength of the projection as indicated by the color bar with black corresponding to absent connections and green for the intrinsic projections that are not included here. The row 
and column ordering was determined by a clustering algorithm based on similarity of the input and output profiles between areas (see the Materials and Methods section). (6) A 
weighted connectivity matrix for the 29 x 29 subgraph. For multiple injections, shading is based on geometric mean values. 



D of nodes of a graph G with node set V is said to be domi- 
nating G, if all elements of V have a link to at least one node 
in D (Kulli and Sigarkanti 1991)- Here, we modify this defi- 
nition slightly by saying that D dominates x% of the nodes of 
G, if an x% of "all" nodes in V are linked to one or more 
nodes in D. The as%= 100% corresponds to "full" domination. 
This definition includes also nodes from D. The minimum 
dominating set (MDS) D min is defined as the one that fully 
dominates G and it has the smallest size (number of nodes). 
For all sets of 2 target area combinations from the 29 target 
areas (406 pairs), 26.6% of them dominate 90-100% of the 91 
areas (see Supplementary Table 5). One pair of areas (81, 7m) 
receives projections from all 91 areas, revealing an MDS 
size of 2. 

A low MDS indicates either a very dense graph or a scale- 
free graph (usually dominated by its hubs) (Barabasi and 
Albert 1999). The in-degree distribution (Fig. 12) is 



"inconsistent" with a scale free graph, as is the fact that 
slightly increasing the size of dominating sets to include 3, 4, 
and more nodes quickly increases their number. For area 
triples there are 69 dominating sets (1.88% of 3654), and for 
sets of 4 areas, there are 1,978 (8.33% of 23 751). Moreover, 
all combinations of 8 sites (out of 29, ~4.29 million) will dom- 
inate at least 90% of all the areas (see Supplementary 
Table 5). As more injections add links but not new nodes, 
they can only enhance these strong domination effects, con- 
firming that the FIN is indeed a dense graph when considered 
in terms of binary connectivity. 

Discussion 

By characterizing the strength of 1,615 identified interareal 
projections, this study provides the most extensive quantitat- 
ive analysis to date on cortico-cortical connectivity in the 
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Figure 12. In-degree distribution. The number of areas projecting to each of the 
target areas of this study. Horizontal dashed line indicates the mean in-degree 57.4. 

macaque. While nearly a thousand of these projections had 
been previously described, published data on interareal con- 
nection weights has been extremely limited. Most of the con- 
nection weights we report are based on only a single 
injection, but we provide evidence that these measurements 
are generally within an order of magnitude of the average 
connection weight for the subregion of the area that was in- 
jected. In addition, we report not only the presence but also 
the weights for 579 NFP. Altogether, this provides a reason- 
ably accurate portrayal of the connectivity profiles for the 29 
injected areas and provides <?29x29 an d ^29x91 subgraphs 
constructed with consistent criteria. Our analysis also pro- 
vides valuable inferences regarding statistical characteristics 
of the entire cortical connectivity graph. To place these find- 
ings into a broader context, several issues warrant discussion. 

Technical Considerations and Limitations 

We used retrograde tracers that show high sensitivity and re- 
stricted uptake zones, the latter making it possible to verify 
their restriction to a single cortical area (for a detailed discus- 
sion of the relative merits of the tracers used here, see Sup- 
plementary material in Markov et al. 2011). While most of the 
injection sites were entirely restricted to the cortical gray 
matter, in a few there was minor encroachment into immedi- 
ately subjacent white matter (see the Materials and Methods 
section and Supplementary Fig. 6). In such cases, a small pro- 
portion of labeled neurons may have arisen from fibers of 
passage. However, this is unlikely to be a major confound, 
especially since the fibers immediately under the gray matter 
mainly arise from the nearby cortex, that is, predominantly 
the same cortical area (Schmahmann and Pandya 2009). 

When combined with close section intervals and full hemi- 
sphere scanning, these tracers revealed many NFP for each 
area injected and enabled estimation of connection weight for 
each projection. However, there are also limitations to the 
methods used. For instance, there are many alternative panel- 
lation schemes for the macaque cortex, some of which are 
finer-grained than the 91-area parcellation we used. Neverthe- 
less, our M132 parcellation has been mapped to the macaque 
F99 atlas, which contains accurate surface maps for 15 other 
published parcellations (Van Essen, Glasser, Dierker, Harwell 
2011), thereby enabling detailed cross-parcellation compari- 
sons. In addition, we report all injection sites by their stereo- 
taxic coordinates in the F99 atlas space. This provides an even 
finer granularity for objectively comparing our connectivity 



data to alternative parcellation schemes and to other data 
modalities such as functional connectivity revealed by func- 
tional magnetic resonance imaging (fMRI) (Vincent et al. 
2007) and structural connectivity revealed by diffusion 
imaging and tractography (Rushworth et al. 2009). It would 
be of interest to determine a quantitative connectivity matrix 
using one or more alternative parcellations applied to each of 
the cortical injection cases. Unfortunately, technical consider- 
ations make this an enormous undertaking that was not feas- 
ible for the present study. 



Sampling Strategies and Statistical Methodologies 

Scannell et al. (2000) emphasized the importance of sampling 
and statistical methodologies in efforts to quantify cortical 
connectivity. Conventional approaches typically involve 
repeat injections into each cortical area and examination of 
connectivity using only a small fraction of the available histo- 
logical sections, and averaging of connectivity data across 
multiple cases. We chose an alternative strategy that includes 
sampling histological sections at high density and across the 
entire hemisphere, thereby capturing as much data as are 
feasible and also reducing variability related to brain mor- 
phology (Falchier et al. 2002; Vezoli et al. 2004; Markov et al. 
2011). We first analyzed the strengths and limitations of the 
single-injection approach by quantitatively analyzing repeat 
injections in 4 areas (14 injections total). We then proceeded 
to analyze connectivity profiles for single injections in each of 
25 areas. 

The repeat-injection analysis indicates that connectivity pat- 
terns across animals are remarkably consistent, except for 
very sparse projections. This was previously shown for central 
visual field injections in areas VI, V2, and V4 (Markov et al. 
2011) and is extended in the present study to area 10, a 
higher order area with about twice as many inputs as the 
early visual areas (Fig. 11). Projections having more than 10 
neurons on average occurred consistently in repeat injections, 
with just 2 exceptions, and therefore can be considered stable 
projections. For very sparse projections, the probability of at 
least one false-negative increases with the number of injec- 
tions (Fig. 5B). False-negatives are also more likely for areas 
such as area 10 that are connected to a larger number of other 
areas (Fig. 5A). The median number of estimated extrinsic 
neurons per hemisphere labeled by one of our repeat injec- 
tions is 187,000 (interquartile range: 98,000-280,000), which 
gives a lower limit to the reliability of the FLNe of <6 x 10 . 

All 4 repeat-injection areas showed a distribution of FLNe 
strengths well fit by a lognormal distribution whose variability 
is described by a negative binomial distribution. The dis- 
persion parameter for each area lies within a limited range 
{9 between 5 and 9). The 95% confidence intervals computed 
for individual injections into areas VI, V2, V4, and 10 include 
97% of the mean values for the repeat injections. These data 
empirically confirm the capacity of our statistical analysis to 
predict the observed variance. The lognormal distribution of 
ordered FLNe values for the 25 single injections is similar to 
those for the 4 repeat-injection areas. Hence, it is reasonable 
to use the negative binomial model derived from the repeat 
injections to estimate the 95% confidence intervals for the 
single injections as well. The lognormal distribution observed 
for all injections indicates that sparse connections are part of 
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an impressively broad and essentially continuous range of 
connection weights. 

Many cortical areas show major internal heterogeneity in 
their connectivity patterns that exceeds the variability 
reported here for repeat injections that were localized to a 
particular areal subregion. This is well documented for the 
peripheral versus central representations of VI, V2, and V4 
(Falchier et al. 2002) and the body representation of F2 
(Johnson et al. 1996; Matelli et al. 1998; Tanne-Gariepy et al. 
2002; Luppino et al. 2003). This cautionary note also applies 
to other topographically organized sensory-motor areas, and 
it might apply to higher-order areas as well (Burman et al. 
2011). This underscores the importance of accurately specify- 
ing the location of the injection sites and the labeled neurons, 
in the context of a consistent database as provided in the 
present study (cf. Fig. 7 and see Supplementary Figs 2 and 4). 

Variability 

There are other sources of variability besides those just dis- 
cussed (statistical fluctuations and regional heterogeneity 
within an area). All 4 of our repeat injections showed con- 
siderable variability, even for projections having moderate or 
strong FLNe values. Some of the variability may be attribu- 
table to imperfect delineation of areal boundaries for the 
various source areas. However, there is essentially no uncer- 
tainty in the boundaries of area VI; yet, the variability of VI 
inputs to V4 exceeds that of many other source areas (see 
Fig. 11C in Markov et al. 2011). Moreover, we saw only a 
modest reduction in variability of the repeat injections in 
areas VI, V2, V4, and 10 (Fig. 30 when the areas were 
grouped into regions (Fig. 3D). 

For the above reasons, we consider it highly likely that 
there are genuine individual differences in connection 
strength that along with the lognormal distribution is a 
characteristic feature of interareal pathways. Whether such 
variability is related to experience during the development 
of the animal and/or individual differences in behavioral 
capacities is an important question that may become accessi- 
ble to future noninvasive imaging studies (see below). We 
also note that the connection strengths reported here are 
(except in one case, see Supplementary Table 2) specifically 
for M. fascicularis; the values might vary for other macaque 
species, including the widely studied rhesus macaque 
(M. mulatto). 

Newly found Projections 

While the NFP are sparser than the known connections, 
overall their weight distribution overlaps with that of the 
known connections (Fig. 9). Further, repeat injections in areas 
VI, V2, V4, and 10 showed that 58% of NFP are present after 
each repeat injection and are therefore classified as consistent. 
We have compared the consistency of the NFP with that of 
the known projections in these repeat injections. This shows 
that at similar weights, the consistency of the NFP is equal to 
that of the known connections (see Supplementary Fig. 5). 

Given the overlap of the FLNe values of known connec- 
tions and NFP, why were the NFP detected here and not in 
previous reports? One factor is that we scrutinized the entire 
hemisphere for labeled cells at high magnification and at 
close section intervals (240 um). Another is that most NFP are 



long-distance connections that may have been missed in 
studies that did not closely scrutinize distant cortical regions. 
A third is that some investigators may have assumed a priori 
that very sparse connections are functionally insignificant 
and could be safely ignored. Finally, as already noted, some 
connections we report as NFP might reflect cross-study differ- 
ences in the designation of areal boundaries in individual 
cases. 

The general finding is that cortical areas containing back 
labeled cells form a relatively continuous field comprising 
multiple areas as illustrated in Figure 7. On the whole, the 
NFP tend to complete the region of labeling formed by the 
known connections. A correlate of this is that in the 29 injec- 
tions, 75% of NFP share borders with areas with known con- 
nections. Hence, in Figure 6 and Supplementary Figure 2, 
there are multiple examples of labeled neurons in NFP that 
could be included in a known connection by shifting a 
border, and in some instances, the shift would result in an un- 
justified increase in the dimensions of the known area. 
However, following random fusing of "all" possible pairs of 
adjacent areas shows that 25% NFP remain (data not shown). 
So while shifting borders could reduce the NFP, it would 
mean that we abandon a consistent atlas with an objective 
parcellation for one where area dimensions are dictated by la- 
beling patterns. Besides decreasing the NFP, such a parcella- 
tion process would also tend to erase another aspect of the 
present findings, namely that there is a drop in the weight of 
connectivity towards the periphery of the region of labeled 
neurons and that neighboring areas share similar connectivity, 
even when labeling levels are very low. 



Reciprocity 

Previous reports provided scant evidence of unidirectional 
connectivity in the cortex (Felleman and Van Essen 1991). 
However, few experimental studies have directly tested uni- 
directionality because it is technically not easy. For example, 
using a bidirectional tracer such as WGA-HRP is not appropri- 
ate. This is because the high levels of recurrent connectivity 
in the source area (Braitenberg and Schuz 1998; Markov et al. 
2011) can lead to labeled axon terminals in the vicinity of the 
retrogradely labeled cells that can be mistaken for long- 
distance anterograde labeling from the injected target area. 
Our evidence for unidirectionality is susceptible to error 
because the reciprocity of connections could only be tested 
using injections in different animals. However, this error 
cannot at present be overcome given that diffusion tensor 
imaging and tractography do not detect direction and existing 
tracing techniques, including BDA and labeled amino acids, 
show bidirectional transport (LeVay and Sherk 1983; da Costa 
and Martin 2011). 

In the present study, reciprocity of cortical connectivity was 
unexpectedly low: There was evidence for unidirectionality 
for 33% of connections, and over 40% of bidirectional connec- 
tions were strongly asymmetrical. However, areal heterogen- 
eity may be a confounding factor for many of the apparent 
unidirectional connections reported in the present study. Posi- 
tive confirmation via superposition of injection sites and retro- 
grade labeling suggests that a minimum of 10% of all cortical 
pathways are genuinely unidirectional. 
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Matrix Density 

The 29 injections reported here revealed many long-distance 
NFP having FLNe values that overlap with the known connec- 
tions and showed similar consistency. This gave rise to a connec- 
tivity density (i.e. the number of existing connections expressed 
as a percentage of the maximum number of possible connec- 
tions) of 66%. Further graph theoretic analysis allowed us to 
infer that the binary FIN is dense and that this is therefore a 
characteristic feature of the macaque cortex. High densities are 
necessarily accompanied by a short average path length across 
the cortex and will impact on large-scale models of cortical 
networks (Young 1993; Sporns et al. 2000; Honey et al. 2007). 

Low cortical matrix densities have been reported in esti- 
mates based on large data sets compiled across studies (e.g. 
Young 1993, who reported a density of ~15%). The lower 
density in these large datasets might suggest that connections 
between specialized cortical systems are uncommon. However, 
subsequent studies have revealed numerous connections 
between functional systems, for example, between early audi- 
tory and visual areas (Falchier et al. 2002). In this respect, we 
found that NFP signficantly span greater distances and more 
frequently link areas in different regions than do the known 
connections (Markov et al. in preparation). If the analysis were 
restricted to previously known connections, the estimated 
density for the entire hemisphere would decrease to 45%, 
similar to that predicted for the visual system alone by 
Felleman and Van Essen (1991) and Jouve et al. 1998. 



et al. 2009; Li et al. 2009). The wide range of densities might 
reflect either differences in thresholding or technical difficul- 
ties (Jones 2010; Jbabdi and Johansen-Berg 2011; Behrens 
and Sporns 2012; Van Essen and Ugurbil 2012). 

While it is widely presumed that functional connectivity 
(i.e. functional correlated fMRI signals across the cortex) are 
constrained by the underlying structural connectivity, the 
nature of this relationship remains to be understood (Bressler 
and Menon 2010). For instance, while there is evidence for a 
good correlation between functional connectivity and strong 
short-distance connections, this is not the case for the more 
numerous long-distance connections (Honey et al. 2009; 
Adachi et al. 2012). A recent modeling study showed that in- 
cluding anatomical directionality and weight led to more rea- 
listic biologically modeled brain dynamics (Knock et al. 
2009). One issue of particular interest is the relative impor- 
tance of projections of different strengths. There are numer- 
ous examples of causative interaction between weakly 
connected areas where the direct cortico-cortical pathway 
may play a functional role (Ekstrom et al. 2008; Wang et al. 
2008). The combination of weight "and" distance of connec- 
tions may be important for understanding large-scale tem- 
poral dynamics (Honey et al. 2007). Sparse long-distance 
connections, conceivably via nonlinear phenomena such as 
"contraction dynamics", may have disproportionate effects 
relative to the far more dense intrinsic and intermediate- 
distance connections (Wang and Slotine 2005). 



Conclusion 

Increasing brain size has been suggested to result in decreased 
connectivity (Ringo 1991). A reduction in connectivity could 
be associated with lower weights of long-range pathways with 
respect to local connectivity or with a reduction in the number 
of pathways (at least in relation to the total number of 
area-to-area combinatorial possibilities). A decrease in weight 
values with increasing size is supported by the comparison of 
interareal connectivity in primate brains of different sizes 
(Palmer and Rosa 2006). Changes in the frequency of long- 
distance connections with changes in brain size would be 
associated with a change in density of the cortical matrix. In 
this respect, it is notable that the mouse cortex has a very high 
density of its interareal graph in comparison to the present 
findings in the larger macaque brain (Wang et al. 2012). 

The human cerebral cortex is 9-fold greater in surface area 
than the macaque and probably contains a maximum of 150- 
200 cortical areas compared with the estimated maximum of 
130-140 in the macaque (Van Essen, Glasser, Dierker, 
Harwell 2011; Van Essen, Glasser, Dierker, Harwell, Coalson 
2011). What are the major differences in cortico-cortical con- 
nectivity in the 2 species? Although the human cortex is not 
amenable to anatomical tracer studies, the emergence of non- 
invasive neuroimaging methods, including tractography using 
diffusion imaging and functional connectivity using resting- 
state fMRI, opens the possibility of learning much about the 
human connectome (Behrens and Sporns 2012; Van Essen 
and Ugurbil 2012; Van Essen et al. 2012). Numerous noninva- 
sive brain imaging studies carried out in primates using brain 
parcellations ranging from 45 to 90 areas report wide ranging 
densities, from 7% to over 50% (He et al. 2007; Chen et al. 
2008; Hagmann et al. 2008; Iturria-Medina et al. 2008; Gong 
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